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1. Introduction 

The vast majority of the existing theoretical and experimental work on spin glasses [[I]] 
has been in regimes where quantum effects can be almost entirely neglected, and both the 
statics and dynamics can be described by classical statistical mechanics models. This 
earlier work has focused on the paramagnetic phase not too far above the spin glass 
transition temperature, T c , and on the phase with spin glass order with temperatures 
T < T c . Thermal fluctuations then dominate, and the typical relaxation frequency, u, 
satisfies Tiu <C ksT, justifying the use of classical models. One important consequence 
of this classical behavior is that the coupling between spin and charge fluctuations is 
rather weak; whether the frozen spins in the spin glass phase are in a metallic, or a 
insulating, host is not of direct importance in theories of the formation of the spin glass 
phase, or of its critical fluctuations. The metallic systems do have a longer range RKKY 
interaction between the spins, but this does not fundamentally differentiate the basic 
theoretical framework used to describe the two cases. 

The situation is quite different when we vary T c 8jS cl function of some tuning 
parameter (e.g. pressure, doping concentration of spins, or disorder) r, and examine 
the regime where T c (r) becomes much smaller than a typical microscopic exchange 
constant |1. In particular, it may be possible to reach a r = r c at which we first have 
T c {r c ) = 0; we will choose an overall shift in the value of r so that r c = 0. For r > 0, 
then, the system is a paramagnet at all temperatures, including at T = 0. Considering 
the situation precisely at T = 0, the ground state undergoes a phase transition, from 
a spin glass (r < 0) to a paramagnet (r > 0), which is driven entirely by quantum 
fluctuations and disorder. Under these circumstances, as is described in more detail 
below, quantum effects play a fundamental role over a substantial portion of the T, r 
plane (and not just at T = 0). Concomitantly, the coupling between spin and charge 
fluctuations becomes much more important. It is now essential to specify whether the 
spin glass and paramagnetic phases are metallic or insulating, as the effective theory of 
spin fluctuations is quite different in the two cases. 

In this paper, we will review recent work on insulating and metallic quantum spin 
glasses, highlighting the physical origin of the differences between the two cases. The 
emphasis will be on the metallic case, primarily because it is of greater experimental 
interest in the context of rare-earth intermetallic compounds described elsewhere in these 
proceedings. A few results described here have not been explicitly obtained before, but 
they are straightforward consequences of work published earlier. 
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2. General considerations 

It is useful to begin by making a few general remarks on the expected phases, and their 
properties, of strongly random electronic systems in three dimensions. For concreteness, 
consider the T = properties of the following Hamiltonian 

H ^ ] tijCi a Ci a ^ ] J^jSi^Sj^ -\- Hi n t, (1) 

i<j,a i<j,fJ- 

where Ci a annihilates an electron on site i with spin a =|, j, and the spin operator 
= c la cr ai3 c ii3 / ^ with o"^ the Pauli matrices. The sites are placed in three 
dimensional space and labeled by i,j, the hopping matrix elements are short-ranged 
and possibly random, and the are random exchange interactions, possibly with 
spin anisotropies. The remainder Hi nt includes other possible short-range interactions 
between the electrons: we will not need to specify them explicitly. 

Figure 1 shows a possible T = phase diagram of H as its couplings are varied. We 
have only identified phases between which it is possible to make a sharp distinction. The 
phases are determined by the behavior of their spin and charge fluctuations. Charge 
transport is characterized by the T = value of the conductivity, a; if o = the 
phase is an insulator, and is metallic otherwise. For the spin sector, the uniform spin 
susceptibility is not a useful diagnostic (as will become clear below), and we distinguish 
phases by whether the ground state has infinite memory of the spin orientation on a site 
or not. The time- averaged moment on a given site is denoted by (Si) (we will generally 
drop vector indices, fj,, v except where needed) and it can either vanish at every site, 
or take a non-zero value which varies randomly from site to site. The average over all 
sites will be denoted by (Si), and will be non-zero only in ferromagnetic phases, which 
we will not consider here. 

The four phases in Figure 1 are: 
(1) METALLIC QUANTUM PARAMAGNET (MQP) 
More simply known as the familiar 'metal', this phase has 

o ^ (Si) = 0. (2) 

The charge transport occurs through quasiparticle excitations (which also carry spin) 
which can be described by a disordered version of Fermi liquid theory ||: the 
quasiparticles still obey a transport equation as in Landau's theory but have wave 
functions which are spatially disordered. However, it is also necessary to include a finite 
concentration of local spin moments to obtain a complete description of the low energy 
spin excitations [|J. These local moments are created by relatively weak fluctuations 
in the disorder, and interact with each other and the conduction electrons via spin 
exchange. Although the spin fluctuations on the local moment sites are relatively slow, 
they eventually lose memory of their orientation, and we always have (Si) = 0. 
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(2) INSULATING QUANTUM PARAMAGNET (IQP) 
This phase has 

a = (Si) = (3) 

and is accessed by a metal-insulator transition from the MQP phase. It is also referred 
to as the 'random singlet' phase ||, the name describing spin singlet bonds between 
pairs of spins. Breaking any of these singlet bonds costs an energy which has a broad 
distribution, with significant weight at low energy. The large density of low energy spin 
excitations means that the spin susceptibility is quite large, and may even diverge as 
T — > (this can also happen in the MQP phase, as argued by Bhatt and Fisher [|J). 
This is the reason we have not used the spin susceptibility as a diagnostic for the phases. 

(3) METALLIC SPIN GLASS (MSG) 
This phase has 

a ± (Si) ± 0, (4) 

and is accessed from the MQP by a spin-freezing transition. The local moments of the 
MQP phase now acquire a definite orientation, and retain memory of this orientation 
for infinite time. The Fermi liquid quasiparticle excitations are still present and are 
responsible for the non-zero a\ the frozen moments appear as random local magnetic 
fields to the itinerant quasiparticles. Alternatively, we may view the spin-freezing 
transition as the onset, from the MQP phase, of a spin density wave with random 
offsets in its phase and orientation, as appears to be the case in recent experiments |J; 
this suggests the name "spin density glass" . Although there are quantitative differences 
between these two points of view, there is no sharp or qualitative distinction and the 
two pictures are expected to be continuously connected. The latter point of view was 
explored by Hertz some time ago, but he did not focus on the vicinity of the T = 
transition between the MSG and MQP phases. An infinite range model with a MSG 
phase was studied in Ref |J with a static ansatz for the order parameter. Here we will 
review the recent complete solution for the MSG phase and the MSG-MQP transition 
in the infinite range model P, 110], the Landau theory for the short range case [IT], EJ 



and the consequences of fluctuations about the mean field theory. 
(4) INSULATING SPIN GLASS (ISG) 
This phase has 

cr = (Si) ^ (5) 

Charge fluctuations are unimportant, and the collective frozen spin configuration is 
expected to be well described by an effective classical spin model. Examples of phases 
of this type may be found in Ref M. 

Let us emphasize that the above discussion was restricted to T = 0, and thus 
describes only the ground state properties of H. Indeed, all of the ground states have 
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been studied earlier, and we shall have little to add to their description here. Our focus 
shall primarily be on the finite T properties in the vicinity of quantum phase transitions 
between the phases. As we will see, this finite T behavior can be quite non-trivial and 
remarkably rich; our results are summarized in Fig |2] and will be discussed in more 
detail later. The possibility that such finite T crossovers may offer an explanation of 
the unusual properties of lightly doped cuprates (and perhaps other strongly correlated 
systems) was discussed in Ref [[TIJ. 

In the following section we will describe solutions of infinite range models of the 
above phases. The insights from these solutions will be used to attack realistic models 
with short-range interactions in subsequent sections. 



3. Infinite range models 



We describe three models, which address different aspects of the phase diagram in Fig 1. 



3.1. Insulating Heisenberg spins 



This subsection reviews the results of Ref |L3 . 

We consider the Hamiltonian H with tij = and H int constructed so that the system 
has exactly one electron per site. The ground state is therefore necessarily insulating, 
and we are considering a model of the IQP and ISG phases. The exchange interactions 
Jij are isotropic in spin space and are independent random variables with = and 
Jfj = Jq/N s , where N s is the total number of sites. As the mean square value of is 
independent of \i — j\, the model has infinite range interactions, and each interaction is 
of order 1 / y/Nl. 

Even this simple model is, however, not solvable. We need to consider a 
generalization of the symmetry of the model from SU(2) to SU(N) (N counts the 
number of spin components of each electron), and also consider the dependence on the 
"length" of the spin at each site (labeled by an integer rib). We then have a family of 
models labeled by the two integers N, rib and the familiar spin-1/2 Heisenberg model 
corresponds to N — 2, rib = 1. It has been shown that the infinite range model is at least 
partially solvable in the limit iV — > oo, rib — > oo but with rib/N fixed at an arbitrary 
value, and also in the limit iV — > oo, rib fixed. 

For N large, the model is in the ISG phase for nb/N > k c , where numerical 
constant whose value is not precisely known. In this phase we have 



\im o (S l (0)-S l (T)}=q EA ^0 (6) 

where r is imaginary time, and q E A the Edwards- Anderson order parameter. The nature 
of the approach of the two-point spin correlator to q E A at large r, and the manner 
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in which q^A vanishes as n^/N approaches k c from above at T = have not been 
determined. 

The IQP phase appears for rib/N < k c , and its properties have been determined a 
little more explicitly. We now have 

(5,(0) ■ S t (r)} ~ - (7) 

T 

for large r at T = 0. The decay of the correlator is quite slow, and as a result, the local 
dynamic linear spin susceptibility, Xl{w) diverges logarithmically as the frequency, u, 
approaches zero: 

~ In ^-A + z|sgn(^) (8) 

In the present model, each spin has the option of forming singlet bonds with a very large 
number of partners, and the diverging local susceptibility of the IQP phase appears to 
be due to strong resonance between different pairings. 

For comparison with the other infinite range models considered below, we note that 
the local spectrum on each site consists of a single doubly degenerate level of a up or a 
down spin, which are then coupled together by the inter site exchange interaction. As 
we will see shortly, this local degeneracy is intimately linked with the diverging local 
susceptibility of the present model. 



3.2. Spins in a metal 

This subsection reviews the results of Refs || |1C|1 . 

We now consider the Hamiltonian H with both Uj and independent, Gaussian, 
random variables, with mean square values independent of The couplings in Hi nt 
are chosen so that the ground state is metallic, and therefore produces a mean field 
theory of the MSG and MQP phases, and of the quantum phase transition between 
them. Unlike the model of Section |3.1| , we have available an essentially exact solution of 
the low energy properties of the present model, including that of the critical properties 
of the MSG-MQP transition, and of the finite T crossovers. 

We will limit our discussion here to T = 0; non-zero T will be discussed later in the 
paper. In the MQP (or disordered Fermi liquid) phase of this model we have 

WP) ■ $(r)> ~ 1 (9) 

for large r. This 1/t 2 decay of spin correlations is characteristic of that found in any 
Fermi liquid. This is perhaps clearer in the Fourier transform to the dynamic local 
susceptibility, which has the linear form 



(10) 
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of the particle-hole continuum. As one approaches the transition to the MSG phase, 
the linear form holds only for very small to, and the dynamic susceptibility is described 
by the crossover function 

X'LH ~ I = (11) 
VA + Vlu 2 + A 2 

The energy A ~ r (recall r is the tuning parameter to the spin glass transition) is 
the crossover scale separating Fermi liquid behavior for w C A, where (|H]) holds, to 
non- Fermi liquid behavior for uo ^> A where x'Li 00 ) ~ sgn(u)-^u. This non- Fermi liquid 
behavior holds at all frequencies at the A = critical point; in imaginary time the 
critical spin correlator behaves as 

1 



<^( ).^(t))~^. (12) 



T 



For r < 0, the system is in the MSG phase. We now have 



c 



(S i (0)-S i (r))^q EA + — (13) 

where c is some constant, and the Edwards Anderson order parameter vanishes as 
Qea ~ r i as one approaches the critical point to the MQP phase. The entire structure 
of the Parisi spin glass order parameter |TJ has also been determined [TT, £|: the results 



closely parallel those of the classical system QXj] with replica symmetry breaking effects 
vanishing as T — > everywhere in the spin glass phase. 

Finally, it is useful to compare the 1/t 2 decay of correlations in the MQP phase with 



the 1/t decay found in Section [D]for the IQP phase. We may interpret this faster decay 
as a consequence of the reduced local degeneracy of each site. Each spin is coupled to 
a metallic fermionic bath, which, loosely speaking, lifts the double degeneracy on each 



site at asymptotically low energies |l4j 



3.3. Paired spins in an insulator 

This is a rather artificial model in the present context, and we discuss it only because 
it is one of the simplest of the infinite range models, and closely related models were 
the first to be solved exactly |15|, 16, |TTJ. Insights gained from their solution were then 



easily transferred H |10[ to the metallic case already described. 

Consider a random Heisenberg magnet with a natural pairing between the spins 
described by the Hamiltonian H p : 

2 

Hp = ^ JpiSn ■ Sj2 + ^ JijabSia ' Sjf,. (14) 

i i<j a, 6=1 

On each site, i, we now have a pair of spins Sn and S^, and this pair interacts with 
an antiferromagnetic exchange J P i > 0. The J P i are independent random variables, 
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but they are constrained to be positive, and have a mean value of order unity. The 
couplings between spins on different sites is however similar to that in H: the Jij a b are 
independent Gaussian random variables with zero mean, and a root mean square value 
of order 1 / ^/W s . 

The key property of this model is in the nature of the local spectrum at each site, 
without including the consequences of the Jijab- The ground state is a non-degenerate 
singlet and it is separated by a gap, J pi , from the excited triplet state. The local 
degeneracy referred to in earlier subsections, has now been completely lifted. It is this 
feature which makes the solution of this model relatively straightforward. A similar 
non-degenerate local ground state, and a gap to the lowest excited state, is also found 
in spin glass models of quantum rotors Jl6| and Ising spins in a transverse field [15]; the 



properties of these models are essentially identical to those described here. 

The local spin correlations in the IQP phase now decay exponentially fast in the 
infinite range model 



(Si(0) ■ 5i(r)> ~ (15) 

where A is an energy scale which vanishes at the transition to the ISG phase. The 
exponential decay is actually an artifact of the infinite range model, and Griffiths effects 
in finite range models lead to a stretched exponential contribution ~ exp(— d \fr\ which 
dominates at long enough times. For the local dynamic susceptibility, these statements 
mean that there is an energy gap A in the IQP phase of the infinite-range model; 
Griffiths effects induce a weak sub-gap absorption of order exp(— c" J\u\) in models with 
finite range interactions. The gap of the infinite range model A ~ r/ln(l/r) as one 
approaches the critical point at r = 0. 

Notice the faster decay of correlations in IQP phase of the present model, as 
compared to those discussed earlier. This is a consequence of the complete lifting of the 
on-site degeneracy. 

Near the critical point, the crossover function analogous to (|TTD is 

Xl{u) ~ sgn{u)Vuj 2 - A 2 6(\uj\ - A) (16) 

So, at the critical point A = 0, we have Xz,( w ) ~ an d spin correlations decay as 1/r 2 . 
In the spin glass phase, the analog of (|TJ) is 



(S i (0)-S i (T))~q EA + - ; (17) 



T 



other properties of the spin glass are similar to those described in Section 3.2. 



4. Order parameter and Landau theory 



This section reviews the results of Refs (1TTII and 
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For the cases of the paired spin model of Section |3.3| , and the spins in a metal 
model of Section it is possible to go beyond the infinite-range model and study the 
corresponding quantum transitions in systems with short range interactions. The basic 
strategy is similar to that followed in the classical spin glass case: introduce an order 
parameter characterizing the important long-time spin correlations, and use insights 
from the solution of the infinite range model, and general symmetry considerations, 
to obtain Landau functionals describing fluctuations in the case with finite-range 
interactions. This procedure was carried out in Ref [TT] for the quantum rotor model 
(whose properties are essentially identical to the paired spin model of Section |3lf ), and 
in Ref || for the metallic spin glass case. No such extension exists yet for the insulating 
Heisenberg spin case (Section [3.1|) , and its development remains an important open 
problem. 

In the following subsections we will (i) introduce the order parameter for quantum 
spin glass to paramagnet transition fill , (ii) obtain Landau functionals for finite- 
range versions of the models of Sections [T3] [[□]] and [T2] ||, (Hi) show that a simple 
minimization of these functionals reproduces the properties of the infinite-range models, 
and describe crossovers for a number of observables. The consequences of fluctuations 
in finite-range models is discussed in Section j|. 



4-1. Order parameter 



We begin by introducing the order parameter for the quantum phase transition [11 



(lb 



Recall that for classical spin glasses in the replica formalism, this is a matrix q 
a,b = 1 . . . n are replica indices and n — > 0. The off-diagonal components of q ab can 
be related to the Edwards-Anderson order parameter, q EA , in a somewhat subtle way 
we won't go into here Jl|, |17 |. In quantum (T = 0) phase transitions, time dependent 



fluctuations of the order parameter must be considered (in "imaginary" Matsubara time 
r), and in the spin glass case it is found that the standard decoupling, analogous to the 
classical case introducing q ab , leads now to a matrix function of two times |18[ which we 
can consider to be 

Q*(x,T lt T 2 ) = £ S?(7i)Sj(T 2 ) (18) 

where A/"(x) is a coarse-graining region in the neighborhood of x. From the set-up of 
the replica formalism it is clear that 

1 



(Si(0) ■ Si(r)) = linj l - E ((Q aa (^n = 0,r 2 = r)» (19) 



n-»0 n 



q EA = lim lim - £ ((Q aa (x, n = 0, r 2 = r)» (20) 

a 
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relating q E & to the replica diagonal components of Q. We have introduced above double 
angular brackets to represent averages taken with the translationally invariant replica 
action (recall that single angular brackets represent thermal/quantum averages for a 
fixed realization of randomness, and overlines represent averages over randomness). 
Notice that the fluctuating field Q is in general a function of two separate times T\ 
and r 2 ; however the expectation value of its replica diagonal components can only be a 
function of the time difference T\ — r 2 . Further, the expectation value of the replica off- 
diagonal components of Q is independent of both n and r 2 fTEfl , and has a structure very 
similar to that of the classical order parameter q ab . One can therefore also obtain q E A 
from the replica off-diagonal components of Q, as noted above for q ab |], l7fl . Let us also 



note for completeness that, unlike the quantum case, the replica diagonal components 
of the classical order parameter q ab are usually constrained to be unity [lj], and contain 
no useful information. 

The order parameter we shall use is Q (x, ti,t 2 ), which is a matrix in a replica space 
and depends on the spatial co-ordinate x and two times r 1; r 2 . However, a little thought 
shows that this function contains too much information. The important degrees of 
freedom, for which one can hope to make general and universal statements, are the long- 
time spin correlations with \t\ — r 2 | ^> r m , where r m is a microscopic time like an inverse 
of a typical exchange constant. As presented, the function Q contains information not 
only on the interesting long-time correlations, but also on the uninteresting time range 
with |n — r 2 | smaller than or of order r m . The correlations in the latter range are surely 
model-dependent and cannot be part of any general Landau action. We shall separate 
out this uninteresting part of Q by performing the shift 

Q ab (x, n, r 2 ) - Q ab (x, n, t 2 ) - C6 ab 5(n - r 2 ) (21) 

where C is a constant, and the delta function 5{t\ — r 2 ) is a schematic for a function 
which decays rapidly to zero on a scale r, m . The value of C will be adjusted so that the 
resulting Q contains only the interesting long time physics: we will see later how this can 
be done in a relatively straightforward manner. The alert reader may recognize some 



similarity between the above procedure, and M.E. Fisher's analysis [19] of the Yang-Lee 
edge problem. In that case, too, the order parameter contains an uninteresting non- 
critical piece which has to be shifted away; we will see below that there many other 
similarities between the Yang-Lee edge and quantum spin glass problems. 

4-2. Action functional 

The action functional can be derived by explicit computations on microscopic models 
or deduced directly from general arguments which have been discussed in some detail in 
Ref IPTf . Apart from a single non-local term present in the metallic case || (see below), 



the remaining important terms are consistent with the general criteria |TTJ that: 
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(i) The action is an integral over space of a local operator which can be expanded in 
gradients of powers of Q evaluated at the same position x. 

(ii) Q is bilocal (i.e. is a matrix) in time, and each time is associated with one of the 
two replica indices (see definition Eqn (jjjf)). These "indices" can appear more than 
once in a term and are summed over freely subject to the following rules before 
summations: 

(a) Each distinct replica index appears an even number of times HJ. 

(b) Repetition of a time "index" corresponds to quantum-mechanical interaction 
of spins, which must be local in time and accordingly can be expanded as terms 
with times set equal plus the same with additional derivatives; it occurs when 
the corresponding replica indices are the same, and only then. 

We now present all the terms, which, a subsequent renormalization-group analysis 
tells us are important near the quantum critical point. This is only a small subset of 
the terms allowed by the above criteria. 

A crucial term is that linear in the order parameter Q. This term encodes the local, 
on-site physics of the spin glass model. To the order we shall consider, this is the only 
term which distinguishes between the finite-range analogs of the models of Section EO 



and For the paired spin model of Section we find the linear term [11] 

d d 



1 ( 

— / d d xdr^2 

Kit J f , 



<9ti dr 2 



_l_ r 



Q m {x,r Xl r 2 ) 



(22) 



The gap in the on-site spectrum for this model tells us that it is valid to expand the 
linear term in derivatives of r, and the above contains the leading order contributions. 
The coupling r will be seen below to be the critical tuning parameter for the transition 
from the spin glass to the paramagnet (it is related to the parameter r, introduced 
earlier, by an additive constant). There is an overall factor of 1/nt in front of this term; 
we have written this factor as a product of two coupling constants, k and t, for technical 
reasons 0| we shall not discuss here. Turning to the metallic model (Section [3. 2D we 



find, instead, the linear term 



— [ d d x\ I dTj2rQ aa (x,T } T) - - f dndr 2 J2 

Kt J [J a TV J a 



g aa (x,r 1; r 2 ) 
T2) 2 



[Tl 



(23) 



Now the coupling of the local spin degree of freedom to the itinerant quasiparticles has 
introduced a long range 1/r 2 interaction in imaginary time: this is clearly related to the 
result @. A term with two time derivatives, like that in (|22~D, is also permitted, but it 
is not as important at long times. 

The remaining terms presented here apply to both the insulating paired spins and 
spins in a metal model. 
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There is a quadratic gradient term 



~' ' a,b 



(24) 



which is responsible for the development of spatial correlations in the spin glass order. 
Clearly, such a term is absent in the infinite-range case, and is special to the more 
realistic short-range case. A quadratic term without gradients 

2 



d d xdr 1 dr 2 ^ [Q*{x, n, r 2 ) 

a.b 



(25) 



is also allowed by the general criteria, but we choose to tune its coefficient to zero by 
using the freedom in (21). As will become clear in the next subsection, this criterion 
is identical to requiring the absence of uninteresting short time behavior in Q. Notice 
again the formal similarity to the theory of the Yang-Lee edge |n|, where setting the 
coefficient of a quadratic term to zero was also responsible for removing the uninteresting 
non-critical part of the order parameter variable. 

Next we consider cubic non-linearities, and the most important among the several 
allowed terms is the one with the maximum number of different time and replica indices: 



-^f d d xdr x dr 2 dn £ Q ab (x, n , r 2 )Q bc (x, t 2 , r 3 )Q ca (x, r 3 , n). (26) 

a,b,c 

This term accounts for non-linearities induced by solely by disorder fluctuations. 

Of the terms with fewer than the maximum allowed number of time indices at a 
given order, the most important one is the one at quadratic order: 
u 



2t 



J dfixdr^u Q aa (x,T,r)Q aa (x,T,r). (27) 



We have ignored vector spin indices (//, v) here, and if these were accounted, we would 
have found two separate terms with the same replica and time integration structure as 
( p7|) JITJ ; however the additional term does not significantly modify the physics, and we 



will ignore it here. The coupling u is the only one responsible for quantum mechanical 
interactions between the spins, and as a consequence, all the time and replica indices in 
(p7|) are the same. 

Finally we have a final quadratic term 



~ J d d x J dr\dr 2 ^ Q aa (x, r 1; Ti)Q bb (x, t 2 , t 2 ), 



. . !2n) 

a,b 



which accounts for the spatial fluctuation in the position of the paramagnet-spin glass 
transition. Recall that the linear coupling f was the control parameter for this transition, 
and a term like is obtained by allowing for Gaussian fluctuations in f, about its 
mean value, from point to point in space. It will turn out that (E3) plays no role in the 
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mean-field analysis in the following subsection. However, it is essential to include ( pp 
for a proper theory of the fluctuations, and its effects are responsible for ensuring that 
the correlation length exponent, u, satisfies the inequality v > 2/d. 

To summarize this subsection, the action functional for the paired spin model with 
short-range interactions is ( P2"|) + fl2~j|) + (|2lp + ( |2"7| ) + (^), while that for the spins in 
a metal case is (|23|) + (|24|) + (|2§ + (|27|) + (|2|). Notice that the two differ only in the 
form of the linear term, which is (^) for the paired spin case, and (^J) for the metallic 
case. 



4-3. Mean field theory 

We will limit our discussion in this subsection to the metallic case. We will review the 
mean field theory for the MQP phase || [Kj and identify the position of its instability 
to the MSG phase. A discussion of the solution within the MSG phase will not be 
presented here, and appears elsewhere ||. 

We Fourier transform from imaginary time to Matsubara frequencies by expressing 
the action in terms of 

-1/T 



dT\dr<iQ 



ab 'x,T 1 ,T 2 )e^ lTl+u)2T2 \ 



(29) 



where we are using units in which fi = ks = 1, and the frequencies, u\, 0J2 are quantized 
in integer multiples of 2nT. Then, we make an ansatz for the mean-field value of Q 
which is x- independent, and dependent only on T\ — T2; within the MQP phase this 
takes the form 



Q ab (x,^ 2 ) = (5 a V+. 2 ,o/T)xL(^i) 



(30) 



where we have used (20) to identify the right hand side as the local dynamic 
susceptibility. Inserting ([3(]) into the action for the metallic case in Section [i~2|, we 
get for the free energy per unit volume T jn (as usual 0, T jn represents the physical 
disorder averaged free energy): 



n t\f 



\uj + r 



K 



K 



-Xl{iu) 



u 

+ 2i 



1 2 



(31) 



Notice that the coupling 1/t appears only as a prefactor in front of the total free energy, 
as the contribution of the 1/t 2 term (^) vanishes in the replica limit n — > 0. The 
value of t will therefore play no role in the mean field theory. We now determine the 
saddle point of (|3l|) with respect to variations in the whole function xl(^), and find 
the solution 



XL{iu) = — \/M + A, 

K 



(32) 
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where the energy scale A is determined by the solution of the equation 



A = r-uT VM + A - ( 33 ) 

Taking the imaginary part of the analytic continuation of ([!2|) to real frequencies, we 
get 

X"M = \= I ( 34 ) 
K v2 a/ A + y/iU 2 + A 2 



which was also the result (|ll|) for the infinite range model. Inserting the solution for xl 
back into (0), and using (|3~3D, we get for the free energy density: 

\ 2" 



T_ 1_ 

n kH 



f E(I-I + a) 3/2 + ^(tEv^^ 



(35) 



The equations (|33|) , fl34|) , (|35|) are key results ||, |10| , from which our mean field 
predictions for physical observables will follow. Despite their apparent simplicity, these 
results contain a great deal of structure, and a fairly careful and non-trivial analysis is 
required to extract the universal information contained within them. 

First, it is easy to note that there is no sensible solution (with A > 0) of ([33|) at 
T = for r < r c where 

r duj f— 2A^/2 

r c = u —yM^— — 36 

where A w is an upper cut-off in frequency. Clearly the system is in the MSG phase for 
T = 0, r < r c , and a separate ansatz for Q is necessary there, as discussed elsewhere 0. 
Let us now define 

r = f — f c , (37) 

so that the quantum critical point is at T = 0, r = 0. In the vicinity of this point, our 
action constitutes a continuum quantum field theory (CQFT) describing the physics of 
the system at all energy scales significantly smaller than A^. The "universal" properties 
of the system are the correlators of this CQFT, and they apply therefore for r, T A u , 
a condition we assume in our analysis below. It is also natural to assume that the 
microscopic coupling u ~ Aj 1 / 2 . We shall, however, make no assumptions on the 
relative magnitudes of r and T. 



Let us now examine the solution of (p3|) under the conditions noted above. Clearly, 
we have from fl3"3] ) and (^?p that 



A + uTVX = r - u \ T V J\u\ + A - / ^J\u\ . (3* 




We have chosen to move the uj = term in the frequency summation from the right 
hand to the left hand side. This permits us to replace the A on the right hand side by its 
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zeroth order result in an expansion in it, as all neglected terms can be shown to be less 
singular as one approaches the critical point. After doing this, we further manipulate 
(BSD into 











f r — 











f du I r — : r — , r 



2\\u 



-u / — = . (39) 
2J\u)\ 



In the first term on the right hand side, the objective has been to always subtract from 
the summation over Matsubara frequencies of any function, the integration of precisely 
the same function; the difference is then strongly convergent in the ultraviolet, and 



such a procedure leads naturally to finite temperature crossover functions pi| . The 
second term also has a subtraction to make it convergent in the ultraviolet. We will 
now manipulate (39) into a form where it is evident that A is analytic as a function of r 
at r = for T / 0. There is a thermodynamic singularity at the quantum critical point 
T = 0, r = 0, but this must disappear at r = for any non-zero T [^jj (see Fig |). We 
use the identity 

^(rb-Tfi)"*^-^ (40) 

to rewrite (39) as 

A + uT\TK = 



+ 




- r^ds (t£ 



du 




s + \uj\ + r J 2ir s + \u\ + r t 
1 1 r 



s + \uj\ f ( s + \u\yj- (41) 

We now evaluate the frequency summation by expressing it in terms of the digamma 
function if>, and perform all frequency integrals exactly. After some elementary 
manipulations (including use of the identity i/;(s + 1) = ip(s) + 1/s), we obtain our 
final result for A, in the form of a solvable quadratic equation for v^A: 



A + uTVA = r(l- + uT 3 ^ (^) , (42) 



where the universal crossover function $(y) is given by 



1 f°° 
Hv) = — / Vsds 

TX A JO 



(43) 
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The above expression for <&(y) is clearly analytic for all y > 0, including y — 0, as 
required from general thermodynamic considerations [21]. Indeed, we can use the above 
result even for y < until we hit the first singularity at y = —2tt, which is associated 
with singularity of the digamma function tp(s) at s — 0. However, this singularity is 
of no physical consequence, as it occurs within the spin glass phase (F ig |2|), where the 
above solution is not valid; as shown below, the transition to the spin glass phase occurs 
for y ~ — mT 1 / 2 which is well above — 2ir. It is useful to have the following limiting 
results, which follow from (H3T) , for our subsequent analysis: 



m = { v 'l/2<(3/2) + 0(y) y _> 

The expression fl34|), combined with the results ( f42|) and (|43| ) completely specify the 
r and T dependence of the dynamic susceptibility in the MQP phase, and allow us to 
obtain the phase diagram shown in Fig 0. The crossovers shown are properties of the 
CQFT characterizing the quantum critical point. We present below explicit results for 
the crossover functions of a number of observables within the mean field theory. A more 
general scaling interpretation will be given in Section |^. 

Before describing the crossovers, we note that the full line in Fig. ||] denotes the 
boundary of the paramagnetic phase at r = r c (T) (or T = T c (r)). This is the only line 
of thermodynamic phase transitions, and its location is determined by the condition 
A = 0, which gives us 

r c (T) = -w$(0)T 3/2 or T c (r) = (-r/w$(0)) 2/3 (45) 

The different regimes in Fig|^ can be divided into two classes determined by whether 
T is "low" or "high" . There are two low T regimes, one for r > 0, and the other for r < 0; 
these regions display properties of the non-critical ground states, which were reviewed 
in Section |2|. More novel is the high T region, where the most important energy scale is 
set by T, and "non-Fermi liquid" effects associated with the critical ground state occur. 
We now describes the regimes in more detail, in turn. 
(I) Low T region above MQP ground state, T < (r/w) 2 / 3 

This is the "Fermi liquid" region, where the leading contribution to A is its T = value 
A(T) ~ A(0) = r. The leading temperature dependent correction to A is however 
different in two subregions. In the lowest T region la, T < r, we have the Fermi liquid 
T 2 power law 

u-nT 2 

ACT) - A(0) = =■ region la. (46) 

6v r 

At higher temperatures, in region lb, r < T < (r/w) 2//3 , we have an anomalous 
temperature dependence 

A(T) - A(0) = w$(0)T 3/2 region lb and II. (47) 
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It is also interesting to consider the properties of region I as a function of observation 
frequency, uj, as sketched in Fig. [3]. At large frequencies, uj ^> r, the local 
dynamic susceptibility behaves like x'L ~ S g n ( u; )\/M') which is the spectrum of critical 
fluctuations; at the T = 0, r = critical point, this spectrum is present at all frequencies. 
At low frequencies, u <C r, there is a crossover (Fig ^ to the characteristic Fermi liquid 
spectrum of local spin fluctuations x'L ~ 00 / V^- 
(II) High T region, T > j\r\/u) 2 / 3 

This is the "non-Fermi liquid" region, where temperature dependent contributions to A 
dominate over those due to the deviation of the coupling r from its critical point, r = 0. 
Therefore thermal effects are dominant, and the system behaves as if its microscopic 



couplings are at those of the critical ground state. The T dependence in ( f47|) continues to 
hold, as we have already noted, with the leading contribution now being A « w$(0)T 3 / 2 . 
As in (I), it is useful to consider properties of this region as a function of uj (Fig |3|). 
For large u {uj ^> uT 3 ^ 2 ) we again have the critical behavior x'L ~ s S n ( Ci; )\/M'; this 
critical behavior is present at large enough uj in all the regions of the phase diagram. 
At small uj {uj ^ wT 3//2 ), thermal fluctuations quench the critical fluctuations, and we 
have relaxational behavior with x'L ~ uj /u l ^ 2 T 3 ^. 
(Ill) Low T region above MSG ground state, T < (—r/u) 2 ^ 3 

Effects due to the formation of a static moment are now paramount. As one approaches 
the spin glass boundary ([4^) from above, the system enters a region of purely classical 
thermal fluctuations, |T — T c (r)| <C -u 2//3 T c 4//3 (r) (shown shaded in Fig |2|) where 

A=(^f) 2 (48> 

Notice that A depends on the square of the distance from the finite T classical phase 
transition line, in contrast to its linear dependence, along T = 0, on the deviation from 
the quantum critical point at r = 0. 

We have now completed a presentation of the mean field predictions for the finite 
T crossovers near the quantum critical point (Fig 0), and for the explicit crossover 
functions for the frequency- dependent local dynamic susceptibility (Fig |3| and Eqns 
(|34|), (|42|), (f43D). We will now consider implications of our results for a number of other 
experimental observables. 



4-3.1. Nuclear relaxation: This was considered in Ref flQfl . The 1/Ti relaxation rate 
of nuclei coupled to the electronic spins by a hyperfine coupling is given by the low 
frequency limit of the local dynamic susceptibility. We have 



A 2 lim 



xUuj) 



A 2 



UJ 



2m/A' 



(49) 
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where A is determined by the hyperfine coupling. The crossover function for \jT\ now 
follows from that for A in (pE2|) , (|43|). In particular, in region II, \jT\ ~ T 1 / 4 



4-3.2. Uniform linear susceptibility: In the presence of a uniform magnetic field, H, 
the action acquires the term M 



| J d d xdr x dr 2 £ Qf v (x,r 1 ,T 2 )H fl H v , 



(50) 



where /i, z/ are spin indices. There are a number of additional terms involving H ||, but 
their linear response is always weaker than that due to (|50|) . Taking the second derivative 
of the free energy with respect to H, we obtain for the uniform linear susceptibility, 

Xu 



11, 10 



Xv 



X6-fVA, 

tK 



(51) 



where Xb is the T, r- independent background contribution of the fermions that have 
been integrated out. The crossover function for Xu now follows, as for 1/Ti, from that 
for A. 



4-3.3. Nonlinear susceptibility: The non-linear response to (ISO) was computed in 
Ref [11]. The nonlinear susceptibility was found to be 

.2 1 



Xnl 



At A' 



(52) 



and its crossover function again follows from that for A. Notice that Xni is proportional 
to the quantum mechanical interaction u, and would vanish in a theory with terms 
associated only with disorder fluctuations. 



4-3.4- Free energy and specific heat: The result for the free energy was given in 
(j35|), and it needs to be evaluated along the lines of the analysis carried out above 
for the crossover function determining A. The specific heat then follows by the usual 
thermodynamic relation. After performing the necessary frequency summations and 
integrations, we obtained the crossover function for the free energy density: 



T{T,r)-T{T = 0,r = 0) 



n 



1 

"^t 



3tt ^VT 
(A-r) 2 



A,„A 2 



2tt 



2u 



4A 5 / 2 
15tt 



where 



2y/2 r°° dVt Q(2y + ^Q' 2 + y 2 ) 



3tt 



y + ^W + y 2 



(53) 



(54) 
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This result for T includes non-singular contributions, smooth in r, which form a 
background to the singular critical contributions. In region II, the most singular term 
is the one proportional to $>p, and yields a specific heat, C v [llOf 



C v C(5/2) r- 

~~rp~ = lb f— VT (55) 



where jb is a background contribution. 

4-3.5. Charge transport: The consequences of the order parameter fluctuations on 
charge transport were explored by Sengupta and Georges |I(J. The quasiparticles are 



assumed to carry both charge and spin, and they scatter off the spin fluctuations via 
an exchange coupling. In the Born approximation, this leads to a contribution to the 
quasi particle relaxation rate, l/r qp 

i r°° ain n 

r qp Jo smh(Sl/T)J A + y / n2 + A 2 



whose T and r dependence follows from that of A. In region II, l/r qp ~ T 3 / 2 . 
5. Fluctuations and scaling analysis 

Fluctuations about the mean field theory just described have been discussed in some 
detail in Refs | ]TT| , |J, and we will be quite brief here. As in Section |4.3| , we will also 
restrict all our discussion to the metallic spin glass case. 

The structure of the corrections can be understood by examining the behavior of 
the action under the renormalization group transformation 

x -> xe~ e t -> re~ zl . (57) 

The mean field theory has dynamic exponent z = 4. 

The system is tuned across the transition by the control parameter r, and by the 
definition of the correlation length exponent, u, it obeys 

dr 1 / _„ x 
— = -r + . . . 58 

at v 

In mean field theory we have v = 1/4. 

The quantum mechanical interaction, u, obeys for weak coupling 

du , . 

Tl = -«.«, (59) 

with 9 U = 2. So the mean field theory in fact corresponds to a fixed point with u* = 0, 
and all u dependent corrections are, in some sense, corrections to scaling. It is essential to 
include such corrections, as none of the remaining couplings involve quantum mechanical 
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interactions between the order parameter modes, and therefore lead to contributions 
insensitive to the value of T. 

The parameter t controls fluctuations of disorder, and for weak-coupling its 
renormalization group equation is 

I = ~ M < 60 > 

with 9 = 2. The value of t was immaterial in the mean- field theory, and it might appear 
from the flow (|60| ) that we can set t = in the analysis of the fluctuations. This not 
strictly correct, as many quantities have a singular dependence on t as t — > 0. As a 
result, t behaves like a "dangerously irrelevant" variable and modifies "hyper scaling" 
relations between exponents. These effects are discussed in more detail in Ref [JTTJ, ||, 
and as they are peripheral to the emphasis here, we will not say more. 
Finally, the cubic non-linearity, k, obeys 

_ = _ K + 9K 3 (61) 

For d > 8, there is a stable fixed point at K = 0, which describes the mean field solution 



discussed in Section [4.3| . However, for the physically interesting d < 8, this fixed point 
is unstable and the system has a runaway flow to strong coupling. So we have no 
systematic way of controlling the corrections due to fluctuations. 

Under these circumstances, the best we can hope to do is to use the insight gained 
from the mean field solution to speculate on some reasonable scaling scenarios. Two 
natural possibilities arise, depending upon whether the fixed point values of u and 
other quantum mechanical interactions are all zero (a static fixed point) or not (a 
dynamic fixed point). The static fixed point will have a 9 U > 0; more generally 9 U is the 
negative of the renormalization group eigenvalue of the least irrelevant of the quantum 
mechanical couplings. Observables sensitive to the order parameter modes will have 
their T dependence controlled by the value of 9 U . Dynamic fixed points have non-zero 
quantum interactions already in the critical theory, and their finite T crossovers were 
reviewed recently in Ref [p2| ; they correspond to the 9 U = case of the finite T behaviors 



of static fixed points, discussed below. 

The remaining discussion summarizes a scaling interpretation of the results of the 
mean-field theory in terms of a static fixed point; we propose that behavior of a generic 
static fixed point will be similar, but with different numerical values for the exponents. 
The crossovers are interpreted in terms of the exponents z, v, 9, 9 U and an order 
parameter anomalous dimension rj. Recall that the mean field theory has z = 4, v = 1/4, 
9 = 2, 9 U = 2, and r\ = 0. Using arguments presented in Ref 0, we can deduce the 
expected behavior of all the phase boundaries in Fig ^ at the inaccessible (but presumed 
static) fixed point, in terms of the above exponents; the results are summarized in the 
caption to Fig Within region II of Fig |2], we obtain the following T dependencies for 
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the observables considered in the mean field theory: 

' rp(d-9-2-2z+r])(l+6 u u)/2z 



•v ^ T( d - - 2 + 1 l)(l+9uv)/2z 

AM 1 

v , ~ r r -(2+z- v -e u )(i+e u u)/z 

Ani 

9jL „ T (d-z-9)/z 
T 

1 ^ r (d-0+2*-2+r,)/2* ^g 2 ) 

The mean field results are obtained by using the mean field exponents and setting d — 8. 

Ref H presented additional speculations on the static fixed point by examining the 
structure of the perturbation theory of the action of Section ^L^. It was found that 
the most singular terms in the perturbation series were identical to those of a much 
simpler classical statistical mechanics problem: that of the Yang-Lee edge in a random 
Ising ferromagnet. This mapping, and some additional assumptions, led to the exponent 
relations 



zv = 1 

1 d-6 + 2-T] 
v ~ 2 



(63) 



As expected, these relations are satisfied by the mean field theory in the upper critical 
dimension d = 8. 



6. Conclusions 

It should be clear that much additional work remains to be done, especially with regard 
to understanding the consequences of disorder-induced fluctuations. Among the issues 
that could be looked at in future work are: 

• Is the quantum critical point for the MSG to MQP transition purely static, i.e. does 
it have no quantum mechanical interactions between the order parameter modes ? 

• A particular feature of our mappings to static critical points is that all exponents 
and crossovers are insensitive to the number of spin components of the order 
parameter, M, at least at all fixed points accessible within perturbation theory The 
mapping to the random Yang-Lee edge, if correct, would also imply an independence 
on M. More recently, Senthil and Majumdar p3"l , have exactly obtained the critical 



properties of a number of random, insulating, quantum spin chains, and found that 
the critical exponents are also independent of M. An important open question, 
then, is whether this independence also holds for realistic metallic and insulating 
quantum spin glasses. 
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• It would be useful to study the classical, random, Yang-Lee edge problem in more 
detail. In particular, although this is a non-trivial problem in d — 2, the critical 
theory is expected to be conformally invariant, and progress may be possible using 
modern methods for such systems. 

• While we have a complete mean field theory for the MSG to MQP transition, there 
is no such theory for the ISG to IQP transition. The IQP to ISG transition is 
not understood even in the infinite-range model, and further progress should be 
possible. 
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METALLIC 
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Figure 1. Schematic diagram of the T = phases of a strongly random electronic 
system in three dimensions described e.g. by the Hamiltonian (^). The average 
moment, (Si), when non-zero, varies randomly from site to site. The conductivity 
is denoted by a. 



25 

A 
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Figure 2. Phase diagram of a metallic spin glass as a function of the ground state 
tuning parameter r and temperature T. In the notation of Fig [l], the T — state 
is a MSG for r < and a MQP for r > 0. The full line is the only thermodynamic 
phase transition, and is at r = r c (T) or T — T c (r). The quantum critical point is at 
r = 0, T = 0, and is described by a continuum quantum field theory (CQFT). The 
dashed lines denote crossovers between different finite T regions of the CQFT: the low 
T regions are la, lb (on the paramagnetic side) and III (on the ordered side), while 
the high T region (II) displays "non-Fermi liquid" behavior. The crossovers on either 
side of II, and the spin glass phase boundary T c (r), all scale as T ~ \ r \zv/(i+0 u v). 
boundary between la and lb obeys T ~ r zv ' . The mean field values of these exponents 
are z = 4, v = 1/4, and 9 U = 2. The shaded region has classical critical fluctuations 
described by theories of the type discussed in Ref M . 
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Figure 3. Crossovers as a function of frequency, u>, in the regions of Fig |2[ The low 
T region is on the paramagnetic side (r > 0). 



